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We study the production of massive fermions in arbitrary vector and axial-vector classical back- 
grounds using effective action techniques. A perturbative calculation shows the different features of 
each field and in particular it is seen that pure temporal axial fields can produce particles whereas 
it is not possible for a pure vector background. We also analyze from a non-perturbative point of 
view a particular configuration with constant electric and axial fields and show that the presence of 
the axial background inhibits the production from the electric field. 
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I. INTRODUCTION 



The production of particles from classical backgrounds has become a very active area of research in the last years. 
We can find it in numerous and disconnected fields of physics such as cosmology Jl],^, heavy-ion collision or even 
£N| \ plasma physics [0 . The pioneer work of Schwinger |i) mainly focused on the production of electron-positron pairs by 
~l • strong electrostatic fields. Since then many other sources of particles creation have been studied in the literature. Thus 
for example we can mention, time dependent gravitational fields Q , varying Planck mass models || , compactification 
of extra dimensions [[|, dilaton field jS], inflaton field, etc. The original Schwinger's work was based on the proper-time 
' technique for the evaluation of the effective action. This method allowed him to obtain an exact non-perturbative 
, result for the total number of particles produced. However, although some other electromagnetic configurations have 
t^J" ' been studied and some particular cases have also been solved || , in general the number of models for which exact 
<^> results can be obtained is very limited. 

In spite of its generalized use, it is probably in cosmology where particle production has been applied more exten- 
sively. Thus for instance, it is believed that the presence of any small anisotropy in the early universe could have 
been erased very effectively by particle emssion processes ||. The theory of reheating after inflation is also based on 
' the resonant production of particles due to the oscillations of the inflaton field (h]] . The reheating of fermions has 
Oh, been considered in ]TT[ | . In addition, the formation of large scale structures in the early universe is closely related to 
the generation from vacuum fluctuation of small seed density inhomogeneities that due to the gravitational instability 
grew to give rise to the current galactic structure. The generation of the vacuum fluctuations can be studied in much 
' the same way as the production of scalar particles in a Robertson- Walker background . Recently the problem of 
particle production has also been extended to the area of the string cosmology in which, apart from the gravitational 
background, there is also an additional scalar field, the dilaton, that can also give rise to the generation of particles 

!& 

In general, the problem of calculating particle production can be approached in two different ways. On one hand, 
the Bogolyubov technique |lj that allows the calculation of the spectrum of the particle produced. It is based on the 
resolution of the harmonic-oscillator equation with variable frequency and only in some particular cases permits the 
derivation of exact results. On the other hand we have the already mentioned effective action technique that allows 
to obtain the total number of particles produced in a much simpler way [|]]lj,[l^] , although the difficulties in finding 
exact results are also present. 

Most of the existing works about the creation of particles from vacuum fluctuations concentrate in the production 
of boson fields in the presence of scalar, vector or gravitational backgrounds. In this paper we will study a different 
source for the production of fermions, it is the presence of general vector and axial- vector backgrounds. We will thus 
extend the Schwinger's work by including the effects of a non- vanishing axial field. Classical axial backgrounds appear 
naturally in modern theories of gravity such as supergravity |l6| or in low-energy string effective actions (lj],|l8j ■ Both 
theories contain torsion (or axion) fields as a fundamental ingredient and, in fact, recently several solutions with 



non- vanishing axial fields have been found in the context of string cosmology [119 20 1. On the other hand, the idea 



of modifying General Relativity by introducing an arbitrary metric connection with torsion is an old one |2l|, and 
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in some sense quite natural from the point of view of the gauge theories of gravity [|22| . This torsion field is coupled 
minimally to fermions by means of its pseudotrace, thus providing a new mechanism for the production of particles. In 
fact in this gravitational theories, torsion would be the dominant mechanism for the production of massless fermions in 
cosmological Robertson- Walker backgrounds. This is due to the fact that when gravity is minimally coupled to masless 
fermions, the theory is conformally invariant. This implies the well known result of absence of particle production. 
However, the presence of additional fields, such as torsion or metric anisotropics breaks that invariance. Our work will 
be based on the effective action method, first from a pcrturbative point of view and then we will study a particular 
case in which a non-perturbative calculation is viable. 

The paper is organized as follows: in section 2 we introduce the lagrangian for the model and also give a brief 
introduction to the effective action technique. In section 3, we perform the perturbative calculation and compare the 
result with the pure vector case. Section 4 is devoted to the particular case with constant electric and axial field and 
a non-perturbative result is obtained in the limit of small axial fields and in section 5 we give the main conclusions 
of the work. Finally we have included an Appendix with some useful formalac of standard perturbation theory in 
quantum mechanics. 



II. MODEL LAGRANGIAN AND THE EFFECTIVE ACTION METHOD 

We will consider the following interaction lagrangian for massive fermions minimally coupled to abelian vector and 
axial-vector fields. For simplicity we will use the left-right notation and at the end we will recover the vector-axial 
fields: 

C = V>(i p - m + ie)ip (1) 

where i fl = i7 p (<9 M + iA^Pl + iB^Pn). As usual the +ie factor is introduced to ensure the convergence of the path 
integral and the left and right projectors are defined as: Pl, = (1 — 7s)/2 and Pr = (1 + 7s)/2. We will use the chiral 
representation for the Dirac matrices in which 75 is diagonal. The coupling constants are included in the own fields. 

Let us now introduce the effective action (EA) for the and fields that is obtained after integrating the 
fermions out: 

(0,t-> oo|0,t-> -00) = Z[A,B] = e lW[A ' B] = N J dipdxpj ^ ^ 

= Ndet(i p - m + ie) (2) 

Here, |0,t — > ±00) denote the initial and final vacuum states that in general will be different due to the presence of 
the external sources. N is a normalization constant that is taken as usual in such a way that 2T[0, 0] = 1, this will 
allow to discard the vacuum divergences as shown below. 

The EA will be in general a complex non-local functional in the external fields. Its real part will contain the 
divergences that will be renormalized by adding suitable local counterterms to the action. The imaginary part will 
be finite and will contain the information about the particle production probabilities. In fact, the probability that 
the vacuum remains stable is given by |(0, t — > oo|0, t — > — 00) | 2 . Therefore the probability that the vacuum decays 
by particle emission will be given by: 

P = l-|(0,t->oo|0,t-» -00) | 2 = 1 - e~ 2 Im W[A,B] _ 2 lmW ^ B ] (3) 

When we only have a vector field such that its corresponding electric field is constant, the probability density per 
unit time and unit volume p can be obtained exactly and the result is given by Q : 

e 2 £ 2 ~ 1 

p = i^E^ e - E ( 4 ) 

71=1 

However for non-constant electromagnetic fields or when the axial part is switched on, the computation becomes very 
involved and it is neccessary to rely on some perturbative method. 

Before concluding this section we will mention that although the notion of torsion appears in strings and supergravity 
in slightly different ways, both can be interpreted as the antisymmetric part of the affine connection in pseudo- 
Riemannian geometry jTsj. Thus, if the components of the metric connection are: r A M „, its antisymmetric part: 

T A M „ = r A M „ — f „„ is known as the torsion tensor. By means of the Einstein equivalence principle, it is now possible 
to minimally couple torsion to fermion fields one gets: 
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C = ^ \d» + %+ l -S^ 5 j i> (5) 
where S p = e^^pT^" is the torsion pseudotrace and fl^, the spin-connection. 

III. PERTURB ATIVE METHOD 

In this section we present the evaluation of the EA in (Q) as an expansion in the external fields, i.e an expansion in 
coupling constants. Let us start by writting: 

W[A, B] = iTr log (i p - m + ie) (6) 

That can formally be expanded as: 

k=l 

where the Dirac propagator is defined as usual by: 

(i - m)-l = [ dqe-^-y) 2 i + ™ . (8) 
" J q z — ra A + ie 

The functional traces Tr are evaluated in dimensional regularization with D — 4 — e and dq — /j, e d ' qj \2ir) D . The 
lowest order contribution in the expansion is given by the two-point terms, i.e: 

W[A,B]W = 1 / d ^ x ^y d p d q^±HL e -^-y)^ Pl+ B V P R ) { + m 2 e-^y- x \fi x P L + B X P R ) (9) 
2 J q z — m z p l ~ m 2 - 

Expanding these terms and defining k = q — p, we obtain the following expression: 

- q^^—^A v y Al + (9Z_*£ + ^-g^AlBZ + (A - B)) (10) 



All the integrals, except for that involved in the term proportional to m 2 , can now be reduced to a common form 
that is evaluated in dimensional regularization: 

h ^-S)tX- m i> = i^m (""^ - D ^ t - £ni - m [ «-/-•) (id 

with I? = rn 2 — k 2 t(l — t). The integral proportional to m 2 is nothing but: 

/ d \q 2 - m 2 )(( q ^k) 2 - m 2 ) - (4^ 2 - ^ 2 > f < 12 > 

The imaginary part of W[.A, B] can be easily extracted from the integrals in the Feynman parameter t. They give 
rise to: 

,i ,i ( 0, to 2 > fc 2 i(l - 1) 

Im / dt(t 2 - t) log(m 2 - k 2 t{\ - t) - ie) = / dt(t 2 - t) I -tt/2, to 2 = k 2 t(l - t) 



-7T, m 2 < k 2 t(l - t) 



7r / 4m 2 (2 Am 2 \ k 2 , s 

zV 1 --w[3 + -w)> ^ >4 (13) 



and in a similar fashion we obtain: 
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Im jf 1 dt{m 2 - k 2 t{l - <)) log(m 2 - k 2 t(l - t) - ie) = -tt-^I- (Jrn 2 - ^\ (14) 

and 

f 1 I 4m 2 /2 A; 2 \ fc 2 

ImJ dtlog(m 2 -k 2 t(l-t)-ie) = -n^l- — {-m 2 - — ), ^>4 (15) 

Putting all the contributions together and changing to the vector and axial-vector fields defined by V = B + A and 
S = B — A respectively, we obtain the final result for the imaginary part in terms of the Fourier transformed fields: 

ImW^lAB] = ^LJ ^ 2 -4m 2 )/l _ ^(-i(l + ^) (F, v (k)F^(-k) 

+ S^(k)S^(-k)) + 2m 2 S^(k)S»(-k)) (16) 

We see that unlike the vector case, the axial contribution to the imaginary part has an additional term proportional 
to m 2 S 2 . This term is prohibited by gauge invariance in the vector case, however as it is well-known it may appear in 
axial theories with massive fcrmions since those theories violate the corresponding gauge invariance. In fact studying 
the divergences that appear in the model we see that they are proportional to the following operators: Fn V F^ v , S^ V S^ V 
and m 2 S^S^ (2^]24|]. The S* 4 operator although having the same dimension does not contribute to the divergent part. 
Therefore it is only neccesary to introduce a kinetic and a mass counterterms for the axial field in order to render the 
theory finite. 

Since the integrand of imaginary part has to be understood as a probability density, it is important to verify that 
it is always positive. As far as k 2 > Am 2 , due to the step function present in ( |l6| ) it is possible to find a reference 
frame in which k = 0. Then we have: 

V) + 2m 2 S ll {k)S»(-k) 

4m 2 )|^| 2 + i(fc 2 + 2m 2 )\V t \ 2 + 2m 2 \S \ 2 > (17) 

From this expression we can extract some of the different features of the production from vector and axial fields. First 
we see that if the axial field is purely spatial in the above reference frame, i.e., So — and we choose Si(k) — Vi(k), 
then the production from pure axial fields is always supressed with respect to the pure vector case. However when 
So(k) = Vo(k) and Si = V\ = then, whereas there is no production in the vector case, it is possible to create particles 
in the axial one. In the massless case both fields give rise to the same amount of particles. 



2m" 



(k S°{k)k S°{-k) - k 2 S^{k)S>*(-k) + S 



k 2 



= o 1 + 2 TT (\Si\ 2 + \Vi\ 2 ) + 2m 2 (\S \ 2 - |5,| 2 ) = -(k 2 - 



IV. CONSTANT ELECTRIC AND AXIAL FIELDS: A NON-PERTURBATIVE RESULT 

In the previous section we have obtained the particle production probabilities up to second order in perturbation 
theory. This is in general a good approximation for small background fileds, however even in those cases, it does 
not contain all the information about the particle production processes. In this section we will study a particular 
configuration of vector and axial fields for which it is possible to find an expression for the imaginary part which is 
non-perturbative in the electric field, in the limit S 2 << E and m 2 S 2 /E 2 < 1. 

Let us start by introducing the operators and acting on states \x) and \p) in the usual form: X^x) = x^\x) and 
Pfj,\p) =Ph\p)' In addition (x\P u \(f)) = id fl (x\(f>). The conumutator is given by [X^^P^] = —ig^ and (p\x) — e lpx / (2n) 2 . 

Following Itzykson and Zuber |26| we recast the Dirac operator in (0) as: 

* p = {F- 4P L - $Pr) (18) 

and taking the transpose we have: 

(P- 4P L - PP H )* = -C(P- A-P R - BP^C- 1 (19) 

where C = ij 2 ^/ is the charge conjugation matrix that satisfies: C-f^C" 1 — —7^ and Cj^C' 1 = 75. The effective 
action in (B) can be written with this notation as: 
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W[A, B] = -»Tr log ({P- 4P L - BP R -m + ie)- 

\ if — m 



m + ie 



(20) 



where we have explicitly introduced the normalization factor N in the last term. The effective action can also be 
written in terms of the transposed operators: 



W[A, B] = -iTrlog (F- 4P L - BP R -m + ief 



If — m + ie 



-iTrlog ((P- JLPr- BP L +m- ie)—± -) 

\ P + m — ie J 



(21) 



adding both expression we get: 



2W[A, B] = -iTrlog (P - J^P L - BP R - m + ie)(P- J^P R - BP L +m- ie) 



P 2 -m 2 + ie 



-iTrlog 
1 



(P„ - A^Pr - B„P L ) 2 - -(A^Pr + B^P L )[^,j v ] + m(4- B)l5 -m 2 +ie 



p2 _ m 2 _|_ i f 

Finally we change to the vector and axial fields, the effective action is then written as: 



2W[V,S] = -iTrlog 



{P» -V^- S^ 5 y - -( V + S^ lb )\r, l v ] - m ^75 -m 2 + iej p2 _ — 2 



m 2 + ie 



(22) 



(23) 



Here V M „ — d^V u — d v V^ and — d^S v — d v S^. Let us take the following background fields: V = (0, 0, 0, Bx 1 ) 
and = (0,0, 0,^) with B and 5* arbitrary constants. This choice correspons to a constant magnetic field B along 
the y axis and a constant axial field S in the z direction. Obvisously the same result will be obtained if we choose the 
fields in different spatial directions, provided they are orthogonal. By means of the Schwinger proper-time integral 
we can write: 

2W[V,S] = -*Tr^°° ^ e -is{m^-U) ^ a ,| e -(P 2 -P 1 2 -P 2 2 -(P 3 -BX 1 -S 75 ) 2 + ^[ 7 1 , 7 3 ]-m5 7 375)| a ,) _ ( X | e -P 2 ( 24 ) 

The action of the traslation operator will simplify this expression: 

(P 3 - BX 1 - S 75 ) 2 = e- 1 ^ {-BX 1 - S^fe 1 ^ (25) 
Let us now introduce complete sets of momentum eigenstates: 



2W[V, S] = -iTr f 
Jo 



ds 



-is(m 2 —ie) 



(2 



nl p is(pl-pl) p ip 3 (p 1 -P 1 )/B e -i(p 1 -P 1 )x 1 



l| p i S (-P 1 2 -(-BX 1 -S 75 ) 2 + iB[ 7 1 , 7 3 ]-mS 7 375)|pl\^ _ [ d P d P r i(p-p)x l n \ r *sP 2 |~ 



(2tt)' 



c (p\e 1 ' 



Performing the integral in the p 3 variable and in the p° and p 2 by means of: 

/CO 
dqe- lsq2 -- 
-oo 

the above expression reduces to: 



2W[V,S] = iTr f 
Jo 



ds 



-is(m —ie) 



B 

8^ 



dp(p\e 



is(-P 2 -(BX + S 7 5) 2 + iB[ 7 1 , 7 3 ]-m5 737 5) 



\p) 



(4tt) 2 S 2 



(26) 



(27) 



(28) 



where for simplicity we have denoted X = X 1 and P = Pi. The integral in p together with the matrix trace can be 
considered as the trace of the evolution operator corresponding to the hamiltonian H = Hq + Hi in ordinary quantum 
mechanics with: 
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, -P 2 - (BX + S) 2 - Ba 2 
H ° ~ 1 -P 2 - {BX - S) 2 - Ba 2 



and 



Hi 



mSa 3 



(29) 



(30) 



with a 1 the corresponding Pauli matrices. The problem of evaluating the effective action is thus reduced to the 
calculation of the spectrum of the H operator. However, since we cannot obtain such spectrum in an exact form, we 
will consider Hi as a small perturbation. With that purpose, we will assume that the contributions to the spectrum 
coming from the perturbation are smaller than the eigenvalues of Ho, i.e, m 2 S 2 < B 2 . We can then apply the standard 
Kato theory for time- independent perturbations in quantum mechanics [ £7J . The Hq operator is made out of two 
shifted harmonic-oscillator hamiltonians with mass M = 1/2, frequency lj = IB and a new coupling to the magnetic 



field. Its spectrum {A^} and eigenfunctions {^"il with {n = 0..oo,i = 1..4} can be easily obtained, they are given 
by: 



,(0) 



i<Pr, 

\ {0} - F + B -7/, (0) - 1 



A (0) - E -B ib {U) - -1 



(0) 



\^=E n + B, if" 



(0) 



3 - V2 





(I) 

i(f>n 

\ k J 



A n,4 = E n - B, C 



(0) 



4 _ s/2 




(31) 



where E n = —2B(n + 1/2) are the energy levels of the ordinary harmonic oscillator and: 



B 



V^2 



"«.! 



H„ 



B 



B 



H n 



B x 



|V. - - 



(32) 



with i/„ the Hermite polynomials. Notice the different functional form of these two functions, it reflects the fact 
that the two harmonic oscillators in (E9|) are displaced in different ways. We will express the spectrum and the 



eigenfuntions of the complete hamiltonian H as perturbative series: \ n ,i — J2 p =o A ^nl an< ^ Vw = J2 p =o ^ni- To first 
order in the perturbation, the spectrum is given by the solutions of the equation Ej : 



dp) 



det KVwI-^lVw) - (A n ,i - Ag)%| = 



(33) 



The solution imply A 
are given by: 



(i) 

n.i 



0, Vn, 2, i.e. there is no first order correction to the energies. The second order contributions 

\(^}\Hi\^)\ 2 



v(2) 



X (0) \ (o) 



(34) 



Although, the Ho spectrum is double degenerated, it is possible to use the above expression valid for non-degenerated 
spectra, since there is no contribution from states in the same multiplet. In any case, the explicit calculation using 
Kato theory yields the same results. The values of the second order perturbation can be evaluated in a straightforward 
way, we get: 



X (2) _ x(2) 
A n,l — A n,3 



2B 



where 



i (2) 
Si,2 



V4 



2B 



K>n = \ {<j>n\<l>n)\ = 



' B 



(35) 



(36) 



with L n the Laguerre polinomials. These polynomials are bounded as n grows, in fact asymptotically we have [f28| : 
L n (x) ~ -l=e x ' 2 x- 1 / i n' 1 / i cos(2Vra-7r/4)+C(n- 3 / 4 ) for x > and therefore rc„ is bounded. In fact for S 2 /B « 1, 



G 



we can expand: = 1 — (S 2 / 'B)(An + 2) + 0(S 4 / B 2 ). Therefore, as expected for small values of the axial field with 
respect to the magnetic one, = 1 is a good approximation. Notice that, although the first correction grows like 
n, the growth must be controlled by higher order terms since the function is bounded. Let us stress that we have 
two different parameters in our problem: on one hand S 2 / B that we assume to be very small and our perturbative 
parameter m 2 S 2 /{2B 2 ). 

Using ( [47l) we can obtain the first correction to the eigenfunctions, we get: 



l<x> 



t in S 
2B 

i/rn S 
2B 



irnS L/,(°)\ 
imS U/,(°)\ 



(37) 



Evaluating now the third order corrections in perturbation theory from (|49| ) they again turn out to be zero. This 
is also the case of the second order corrections to the eigenfuntions, i.e. from (^) we get \ipnl) = Vn, i. Thus 
we can calculate fourth order corrections to the energies from ( |49| ) which are non-vanishing. Finally, we see that 
the fifth order correction again vanishes. In conclusion, our results for the perturbed spectrum, up to sixth order in 
perturbations is given by: 



A„,i = A„, 3 = B ( -2 ( n 
A„,2 = A„. 4 = B ( -2 I n 



f 



m 2 S 2 
~2B 2 ~ 
m 2 S 2 
~2B 2 ~ 



m 4 S 4 
8B 4 
m 4 S 4 
~~8B r 



O 



O 



3 s 6 



B e 



">S 6 



B 6 



(38) 



Once we know the perturbed spectrum we can readily calculate the traces in (|2q). We will only perform the Dirac 
traces, but not the functional trace that is equivalent to the integration J d x. Thus we obtain the result for the 
effective lagrangian w that as expected does not depend on x: 



2w[V,S] =4i f 
Jo 



ds 



—is(r 



<0 



B 
8^2 



cos \sB 1 



m 2 S 2 m 4 S 4 



2B 2 



8B 4 



O 



>S 6 



B 6 



,t«(-2B(n+l/2)) 



{Att) 2 s 

Finally performing explicitly the addition of the series in n, we obtain 

sB(l + ^r - ' 



(39) 



2w[V, S] 



ds 



-is{m —ie) 



COS I 



B- 



4 

— 



o 



(*))) 



sin(s_B) 



(40) 



It can be seen that the result is purely real |29J, i.e, there is no particle production in the presence of constant magnetic 
and axial fields. In addition, the integral in s is divergent in the ultraviolet limit s — > 0. As usual these divergencies 
has to be removed by adding suitable counterterms. In order to obtain them, let us expand the part of the integrand 
in brackets around s = 0, we have: 



cos 



B- 



1 



1W 



111 S 
8B 4 



o 



(*))) 



sin(si?) 



B 



1 

~Bs 



sB 
~2~ 



sm 2 S 2 
2B 



O 



m 6 S 6 



+ C(s 3 ) (41) 



The 0(s 3 ) give rise to finite contributions when integrated. The first term is exactly cancelled by the 1/s substraction. 
But we will have to include new counterterms proportional to s(B 2 + m 2 S 2 ). Notice that these are exactly the same 
operators that we found in the perturbative calculation in section 3. In fact there is no contribution from S 4 operators 
due to an exact cancellation of the different cuartic terms in the expansion in (|4l]). Since the theory is renormalizable, 
it is neccessary that the same kind of cancellation operates on the higher order terms in the r.h.s. of the above 
expression. 

The presence of the axial field is known that does not introduce any gauge anomaly in the electromagnetic current 
p0| ], therefore the effective action will be gauge invariant. As a consequence it should be built out of scalar and 
gauge invariant functions. In our case, with constant magnetic and torsion fields, the only possibilities are: F^F^ = 
2{B 2 — E 2 ) and F^ V F^ V = (4B ■ E) 2 . Since the second term vanishes in our case with constant magnetic or electric 
fields, the effective action is invariant under the transformation B —> —iE |29| |. In this way we can convert our results 
for constant magnetic field to pure constant electric fields. Taking into account that now sin(s£>) — > —ismh(sE) 
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and cos( Sj B(1 + m 2 S 2 /(2B 2 ) - m i S 4 /(8B i ))) -> cosh(s.E(l - m 2 S 2 /(2E 2 ) - m 4 S 4 /(8E A ))) and integrating using the 
residues technique we can obtain the expression for the imaginary part of the effective lagrangian: 



f °° (— l) n 
p = 2lmw\E,S] = —^E 2 V V I cos | utt 
47r 3 n 2 

n=l 



m 2 S 2 
~2E 2 



1 



m 4 S 4 
8£ 4 



O 



m 6 S 6 
~E~s~ 



(42) 



When S = we recover the original Schwinger result in (Q) . Notice that in the absence of electric field there is no 
particle production even with non-vanishing axial field. When the mass is zero, we recover the usual result, i.e, axial 
fields only contribute in the massive case. In Figure. 1 the upper curve represent the probability density p as a 
function of S for E = 0.1m 2 , this value for the electric field ensures that the condition S 2 < E is satisfied for all the 
values of S in the plot. 



2e-18 




0.06 



Figure 1. -Probability densities p in units m 4 versus axial field S in units m, for an electric field E = 0.1m 2 . The upper 
curve represents the perturbative calculation up to sixth order. The lower curve is the resummation of the perturbative 
series estimated in the text. 



We can try to extend further the above result by means of the following observation. Since the only possible 
divergences in our model are those mentioned before, it is neccessary that, when expanding the cos function in ([IT 
the higher order terms in S cancel, this implies: 



sB 2 



m 2 S 2 m 4 S 4 



2B 2 



8B 4 



1 q2\ 



s(B 2 + m 2 S 



(43) 



Therefore we can obtain the complete result for the effective action to all orders in perturbation theory in the limit 
S 2 << B. Performing the rotation to electric fields, the result is given by : 



p = 2lmw[E,S] = -^E 

47T J 



2 E 

n=l 



(-1)' 



COS 717T 



1 - 



n 2 S 2 
~EP~ 



(44) 



In Figure. 1, this probability is represented by the lower curve. 

We see that the effect of the axial field is to suppress particle production. Eventually it could make it to vanish. The 
point of vanishing p should indicate the breakdown of the perturbative approximation since the probabilities should 
be positive. In order to check whether the perturbative calculation is valid up to the point in which P vansihes, we 
will study the convergence of the perturbative series. In particular, there is a general result due to Kato M that 
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states that if there are two non-negative constants a and b such that: ||-ffi|| < a||Ho|V')ll + b\\ IV 7 ) 1 1 f° r ai l IV') € D(Hq) 
and if the operator H± is bounded then the eigenvalues perturbative series is absolutely convergent if 2||i?i|| < d. 
With d being the distance from A" i to the rest of the spectrum of Ho- In our case, we take a = and provided 
{\ipn, »)} is a complete set of eigenfunction we can expand Hi\tp) = Hi J2 n » Cn,i|V , n,<)- It is then easy to show that: 
H-^ilV')!! 2 = m 2 S 2 J2 n i C 2 ^ = m 2 S 2 . Therefore if we take b = mS, since Hi is bounded (it is constant) and in our 
case d = 2B, we have that the series is absolutely convergent for m 2 S 2 / (2B 2 ) < 1/2. Transforming to the electrostatic 
case, if we look at the plot we can realize that for those particular values, the point in which the curve crosses the 
axis is in fact approximately signalling the breakdown of the perturbative series. 



V. CONCLUSIONS 



In this work we have studied the production of massive fermions from a classical vector and axial- vector background. 
Using a perturbative evaluation of the effective action we have obtained the contributions to the imaginary part of 
the effective action up to second order in the external fields. We have shown that in the reference frame for which 
k = 0, the production from purely spatial axial fields is suppressed with respect to the that of the vector background. 
In addition, for purely temporal axial fields it is possible to create particles whereas this is not the case for vector 
fields. 

In the particular case of a small constant axial field and a constant electric field, it is shown that a non-perturbative 
calculation can be carried out when those fields are orthogonal. In this case, it is shown that in the massless limit 
the axial field does not affect the production from the electric field. However, in the massive case, the presence of the 
axial background inhibits such production. 

Finally, let us compare these result with the anisotropy damping phenomenon at the Planck era. As is well-known, 
the presence of small anisotropics in the early universe can be damped in a few Planck times due to the backreaction 
of the particles produced on the geometry. An interesting possibility is that a similar mechanism could take place in 
the presence of some primordial torsion field. In this case, since torsion can be generated by the initrinsic spin, this 
could happen when the fermions are produced in some configuration such that the total spin angular momentum of 
the system did not vanish. The use of effective action methods for fermions could be extended to the production of 
higher spin fields, such as gravitinos in a straightforward way. In addition, it is also interesting to study not only the 
particle production rate derived from the effective action, but also the spectra and angular distribution of fermions 
produced. This could be approached by means of the traditional Bogolyubov technique. Work is in progress in this 
direction |3^| . 



VI. APPENDIX 



In this Appendix we summarize the main formulae of the standard perturbation theory used in the text. Let us 
assume that the hamiltonian of theory can be decomposed in: 

H = H + H 1 (45) 

As shown in the text, we denote by \^ \ and the eigenvalues and eigenfunctions of Hq, which are assummed 

to be known. The unperturbed spectrum is assumed to be non-degenerated and discrete. The eigenfunctions and 
eigenvalues of the complete hamiltonian H are expanded in series: \ n ,i = J2 P =a ^n,i an< ^ iVV-t) = S P =o \^n,i) wnere: 

and 

IV© = S n ,i (hi\^) -f>g|V>i7 fc) >) (47) 

with 

S ^ ~ 2. A (o) _ A ( ) ( 48 ) 
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Some simplified formulae for the lowest order terms in the spectrum are given by: 



y( 2)_ v i(V'Sigi^S)r 

™' 1 2^ ,(0) _ A (o) 
\(3) _ /,/,(!) I rr 1 Wu(l)\ 

xS = <<]|^-a£M>-aS<<1i<]) 

AS = Cl^i - ASl«) - 2Ai>e(<^S) - AS(^]|^)) (49) 
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